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We have investigated the effects of long-range Coulomb interaction on the topological super¬ 
conducting phase in a quasi-one dimensional semiconductor wire, proximity coupled to a s-wave 
superconductor using the exact diagonalization approach. We find that in accordance with previous 
studies the addition of Coulomb interaction results in an enlargement of the region of parameter 
values where topological superconductivity can be observed. However, we also find that although 
the interaction decreases the bulk gap for values of the magnetic field close to the phase transition 
point, for moderate magnetic fields away from the transition point, the interaction actually enhances 
the bulk gap which can be important for observation of topological superconductivity in this system. 


Majorana fermions (MFs) [I| have gained considerable 
attention in the past few years due to several proposals 
for the existence of Majorana modes in semiconductor 
systems [H as an elementary excitation. There were 
several experimental attempts to observe the MFs in 
semiconductor systems |1M14||. One of the most promis¬ 
ing candidate for realization of the MFs ii is the ob¬ 
servation of the topological superconducting phase in a 
one-dimensional (ID) semiconductor wire proximity cou¬ 
pled to a s-wave superconductor and with large Rashba 
spin-orbit (SO) coupling [l^. Considering a semicon¬ 
ductor system with high Lande g factor and applying a 
magnetic field perpendicular to the Rashba SO coupling 
direction, the Kramers degeneracy can be lifted by induc¬ 
ing a large gap between the two states. By tuning the 
chemical potential of the system in the gap region, the 
system effectively becomes spinless and supports MFs at 
the edges of the wire similar to Kitaev’s p-wave supercon¬ 
ductor chain model [l^. It should be noted that the MFs 
can also be found in two-dimensional (2D) quantum wires 
where several subbands (channels) are occupied [iMl, 
provided that the number of occupied subbands is odd 
and the width of the wire is smaller than the supercon¬ 
ducting coherence length. Tilting of the magnetic field 
and the issue of orbital effects in 2D semiconductor wire 
have also been addressed [2^. 

The major attraction for finding the MFs is due to 
their potential use in topological quantum computation 
1^ . The non-Abelian nature of the exchange statis¬ 
tics of the MFs makes fault-tolerant topological quantum 
computation feasible, which avoids the issue of decoher¬ 
ence in such a system. While for a strict ID wire the 
exchange statistics is ill-defined, one can form networks 
of ID quantum wires and then move and exchange 
the MFs using closely spaced electronic gates. Therefore 
the quasi ID wire systems where the MFs can be ob¬ 
served, has both theoretical and technological interest. 

Most of the recent studies have focused on the char¬ 
acteristics of topological superconductivity and the exis¬ 
tence of the MFs in semiconductor wires without the in¬ 
teraction between the electrons in the wire. There were 
a few attempts to include electron-electron interaction 
in a strictly ID wire using bosonization, density matrix 
renormalization and Hartree-Fock methods [l^, and 


also in a 2D wire using the mean-field approach [2^. The 
effect of the electron-electron interaction on charge jumps 
(which can be used to observe the MFs) in a ID wire 
has also been discussed It was observed that the 
interaction renormalizes the parameters of the system, 
namely, it enlarges the region of the chemical potential 
and the magnetic field strength, where topological super¬ 
conductivity can be observed. For the ID wire it was also 
found that interaction suppresses the bulk gap, and can 
eventually destroy it for high values of the interaction 
strength. However, the effect of the long-range Coulomb 
interaction on topological superconductivity has not been 
addressed as yet. 

In this work we analyze the effect of long-range 
Coulomb interaction on the topological superconducting 
phase in a quasi ID semiconductor wire, proximity cou¬ 
pled to a s-wave superconductor. In order to do that we 
employ the exact diagonalization scheme to obtain the 
energy eigenstates and the wave functions for the system 
with non-constant electron numbers with odd or even 
parity. We find that in accordance with previous stud¬ 
ies, the addition of the Coulomb interaction results in 
an enlargement of the region of parameter values where 
topological superconductivity can be observed. We also 
find that although the interaction decreases the bulk gap 
for values of the magnetic field close to the phase tran¬ 
sition point, for moderate magnetic fields away from the 
transition point it actually enhances the bulk gap, which 
can be important for observation of the topological su¬ 
perconductivity in this system. The exact diagonaliza¬ 
tion procedure employed here is quite general and can 
be used to address the issue of orbital effects and multi¬ 
channel 2D quantum wire for observation of topological 
superconductivity in the presence of the Coulomb inter¬ 
action between electrons. 

We consider a 2D semiconductor wire with hard wall 
confinement and a strong Rashba SO coupling M in an 
applied magnetic field. As has been considered before, 
the semiconductor ID wire is proximity coupled to a s- 
wave superconductor. We take the wire to be situated in 
the xy plane, with ^ Ly, where L,,, and Ly are wire 
sizes in the x and y directions respectively. Without the 
superconducting pairing potential the Hamiltonian of the 
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system is 


^ = E^SP + 2 E^*r ( 1 ) 

i i^j 

Here = e^/e |r, — rj is the Coulomb interaction term 
and Tigp is the single-particle Hamiltonian, which can be 
written in the form 


Hsp — -1- 5go + Hz, 

( 2 ) 



•w - 2m ^ + ^c(2^-y). 

( 3 ) 

Hso = [n X cr]^ , 

( 4 ) 

Hz = • cr. 

( 5 ) 


TLt^ is the kinetic energy plus the confinement poten¬ 
tial for the system with the chemical potential where 
n = p-b (e/c)A is the canonical momentum, A = —B^y 
is the vector potential, m is the effective mass. The con¬ 
finement potential is V{x,y) = 0, when —L^I2 < x < 
L^I2 and —Lyl2 < y < Byl2, and V{x,y) = oo oth¬ 
erwise. BgQ is the Rashba SO interaction term, with 
the SO coupling strength a. The Rashba SO coupling is 
considered to be present due to the confinement or the 
external electric field which create an asymmetry in the 2 
direction. Finally, i® Zeeman energy term, where 
g is the Lande g factor for the semiconductor. The mag¬ 
netic field is taken to lie in the xz plane with components 
B = (H sin 6,0, B cos 6). By taking as the basis states the 
eigenstates of the when 5 = 0 , namely 



ix + L^/2)-K 


n-y {y + Ly/2)'n 


( 6 ) 

( 7 ) 


the Hamiltonian m can be cast into the second quantized 
form 


E i'^^\'^Sp\^'s')cisCu's' + 

nsn' s' 

\ E {niSin2S2\V^2\n'is'in'2s'2)y. 
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( 8 ) 


where the pairing potential strength A is taken to be 
real. 

Instead of considering the Coulomb interaction at the 
mean field level [ 2 ^ where one writes the complete Hamil¬ 
tonian 'HpgQ = % + "Hgp in the Bogoliubov-de Gennes 
form, we directly treat the Coulomb interaction. In order 
to do that we use the exact diagonalization procedure 
to diagonalize Bpsc even and odd sectors for small 
system sizes. For example, for the odd sector we diag¬ 
onalize 'Hpsc ^ system with non-constant number of 
electrons, namely 1,3, ...Ag electron number basis. A 
similar procedure is employed for the even sector as well. 
This gives us the possibility to obtain the low-lying en¬ 
ergy states and the wave functions both for even and 
odd sector. This procedure can be used to investigate 
both the existence of the topological superconductivity 
and the effect of interaction on the pairing induced bulk 
gap. We use two complementary approaches to identify 
the topological superconducting phase described previ¬ 
ously [24| . If the electron number parity is a conserved 
quantity (which is the case for the Cooper pairing po¬ 
tential), then it is obvious even from the Kitaev model 
that the two degenerate ground states of the topological 
superconducting phase have different parity. Therefore 
the first notion of the topological superconductivity can 
be obtained by considering the value of 

/S.E = , (10) 

where and are the ground states in the odd 

and even sector respectively. In ordinary superconduc¬ 
tors the ground state is unique with integer number of 
Cooper pairs and therefore it has even parity. Hence A5 
is finite in ordinary superconductors. In the topological 
superconducting phase, when two MFs are exponentially 
localized at the two ends of the wire the AA is expo¬ 
nentially small. Therefore, A5 serves as some kind of 
order parameter to distinguish between the topologically 
trivial and non-trivial phases. 

In the second approach we calculate the Majorana 
wave functions from the obtained even and odd parity 
ground states. Let |0) and |1) be the even and odd par¬ 
ity ground states respectively. As was shown in the Ki¬ 
taev model the odd sector state | 1 ) is obtained by adding 
one non-local fermion, which is composed of two MFs lo¬ 
calized at the two ends of the wire, to the even sector 
ground state | 0 ), i.e., | 1 ) = /^| 0 ), where the fermion op¬ 
erator / = 5 ( 71 - 1 - * 72 ) and 7 i, 72 are the MF operators. 
Here the MF operators satisfy the well known relations 
7 ^ = 1 and { 7 q, 7 ;,} = Using these operators we 

can write 


where for brevity n = {n^,, n^} has been introduced and 
s denotes the spin quantum number of the particle. The 
proximity induced superconductivity potential can be 
written directly in the basis introduced above and has 
the form 


|l)=7i|0) = -*72|0). ( 11 ) 

Using the creation and annihilation operators c|js and 
for the basis dH)-©, we can expand the MF operators 
7 i ,72 in the form 


n 


la 


E + 



( 12 ) 
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where are the expansion coefficients. By noting that 
{cL>7a} = and using ((TTl) we get 

= (0 |cL| 1) + (1 |cL| 0), (13) 

‘Pn} = -*(0|cL|l)+*(l|cL|0). (14) 


After obtaining numerically these expansion coefficients 
the probability distribution of the MFs can be obtained 

via the relation p(“)(x,y) = X)* E„ (a:)<(>„. iv) 


We now discuss our results for small number of elec¬ 
trons in a semiconductor wire proximity coupled to a 
s-wave superconductor using the exact diagonalization 
technique. The calculations were performed for the InAs 
semiconductor wire with following parameters: m = 
0.042mQ, where TOq is the bare electron mass, g = —14, 
e = 14.6 [ 2 ^. The SO coupling strength was taken in 
all calculations to be a = 45 meV • nm and the supercon¬ 
ducting pairing potential strength A = 0.225 meV. The 
wire sizes in all our calculations are = 3000 nm and 
Ly = 150 nm. For the even sector we have considered up 
to eight electrons and seven electrons for the odd sector. 
This means that the number of electrons for the even 
sector in the many-body basis takes even values in the 
range from 0 to 8. Similarly for the odd sector it takes 
odd values in the range from 1 to 7. In order to achieve 
numerical convergence in our exact diagonalization calcu¬ 
lations we have introduced a gap between the basis states 
<7 and >7. In a real system, this kind of gap can 
be achieved by adding a periodic potential to the system. 
The addition of the periodic potential can be regarded as 
a spatially varying chemical potential, and therefore it is 
desired that the periodic potential strength be smaller 
than the range of the chemical potential where topolog¬ 
ical superconductivity can be observed. Comparing this 
situation to the lattice model considered previoulsy 
it is natural to take < N^, where is the maximum 
number of electrons in the system, so that every site is at 
least partially filled, which is the case in the density ma¬ 
trix renormalization group studies [ 2 ^. In this work we 
only discuss the 0 = 7r/2 case, which means that the mag¬ 
netic field is aligned along the wire axis, and therefore we 
do not consider the orbital effects here. We also consider 
only the first transverse mode, i.e., we take = 1 in 
all our calculation and therefore do not consider multi¬ 
channel effects. These two issues will be addressed in our 
future studies. 

In Fig. [T] (a) the dependence of the absolute value 
of the difference between the ground state energies in 
odd and even sector {AE defined in (fTU]) l on the mag¬ 
netic field B is shown for the chemical potential p = 0.4 
meV. Both cases of with and without the Coulomb in¬ 
teraction are presented. The value of ^ = 0.4 meV 
corresponds to the energy of the first transverse mode 
Hy = 1. For the quasi ID wire the transition between the 
topologically trivial and non-trivial phases occurs when 
14 = J(e„^ - m)^ + A2 [13, where I 4 = gg.^B/2 and 


= h‘^Tr^ny/2mLy is the transverse mode energy. From 



FIG. 1: (a) The dependence of the absolute difference between 
the energies of the ground states of odd and even sector on the 
magnetic field B. The chemical potential is taken to be p = 
0.4 meV, which corresponds to the transverse mode energy 
for Uy = 1. (b) Same as in (a) but for the dependence on the 
chemical potential for B = 0.7 T. (c) The dependence of the 
bulk gap value (difference between the first excited and the 
ground state energy) for even and odd sector on the magnetic 
field B for the chemical potential g = 0.4 meV. (d) Same as 
in (c) but for the dependence on the chemical potential for 
B = 0.7 T. 


Fig. m (a) it is clearly seen that with an increase of the 
magnetic field strength, AE decreases for small values of 
the magnetic field, becomes zero at some point and then 
starts to increase. The zero point corresponds to a phase 
transition, where below that point the system is in the 
topologically trivial phase, whereas above that value it 
is in a topologically non-trivial phase. Due to the finite 
size effects AE is not exponentially small in the topo¬ 
logical superconducting phase. For the non-interacting 
case the transition occurs at B^ = 0.64 T which is close 
to the value of = 0.56 T for the ideal system calcu¬ 
lated according to the relation above. Clearly, inclusion 
of the interaction results in the lowering of the value of 
which means that it extends the region of the mag¬ 
netic field values where the topological superconductivity 
can be observed [13, [13 • 

In Fig. [ 1 ] (b) the dependence of AE on the chemical 
potential g is shown for B = 0.7 T. As can be seen here, 
AE becomes zero at two points, as expected, because 
for an idealized system the topological superconductivity 
should be observed in the range e„ — < g < 

e„ -I- — A^. For an ideal system and for the values of 

parameters considered here we get the range 0.23 meV < 
g < 0.57 meV, which is quite close to the values observed 
in Fig.[T](b). Similar to the case of Fig. [T] (a) and also to 
the previous studies [13, HI], we observe from Fig. [T] (b) 
that inclusion of the interaction results in a broadening of 
the range of g where topological superconductivity can be 
observed, except that in the present case the broadening 
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is in the lower values of the chemical potential. 

In Fig. [T] (c,d) the dependence of the energy difference 
between the first excited state and the ground state for 
both even and odd sectors on the magnetic field B and 
the chemical potential ^ is shown. As can be seen from 
these figures, in the topological superconducting phase 
the gap in the even sector is smaller than in the odd 
sector. Therefore this smaller gap corresponds to the 
bulk gap of the topological superconducting phase. In 
Fig. [T] (c) this bulk gap decreases with increasing mag¬ 
netic field for the non-interacting case. This is due to the 
fact that increasing the magnetic field aligns the spins of 
the electrons opposite to its direction, therefore reducing 
the pairing potential and destroying the superconducting 
phase. Inclusion of the interaction has two kinds of effects 
on the bulk gap. For magnetic fields close to the transi¬ 
tion point value the bulk gap is reduced. Therefore 
although the interaction broadens the region of the su¬ 
perconducting phase, it lowers the bulk gap, which makes 
the observation of the topological superconducting phase 
troublesome [2^. Increasing the magnetic field strength 
further results in the dependence of bulk gap to become 
almost flat in case of the interaction. In Fig. [T] (c) for 
B = IT the bulk gap is already larger for the interacting 
case than that for the non-interacting case. As was shown 
previously [l^, for even higher magnetic fields the in¬ 
teraction eventually suppresses the superconducting bulk 
gap. Therefore at the topological superconducting phase 
there is some middle region of the magnetic field values 
where the interaction favors the superconducting phase 
and enhances the bulk gap value. This finding can be 
important for possible observation of topological super¬ 
conductivity in this system. Finally, in Fig. [T] (d) we see 
that for B = 0.7 T, the interaction lowers the bulk gap 
for all values of where the topological superconductiv¬ 
ity can be observed. 

In Fig. [5] (a,b) the difference between the single-particle 
densities of odd and even sector many-body state is 
shown without and with the inclusion of the Coulomb 
interaction. For all figures in Fig. [2] the parameter val¬ 
ues /r = 0.4 meV and B — 0.7 T (I 4 = 0.28 meV) have 
been used. In the ideal system the single-particle den¬ 
sities should be the same for odd and even sector [^ . 
and therefore this density difference should be zero. In 
our calculation for the cumulative particle number dif¬ 
ference through all the points of the wire and for the 
non-interacting case we get dN = 0.1 particle difference 
which is also quite close to the value obtained previously 
for the ID wire [l^. When the Coulomb interaction is 
includede, this number decreases to the value 6N = 0.03. 
While it was predicted that for the non-interacting case 
the charge due to the MFs will have an oscillatory behav¬ 
ior in the wire and will be spread uniformly, as is shown in 
Fig. [5] (a) the difference between the single-particle den¬ 
sities of the odd and even sector is not oscillatory and 
is mostly localized at the center of the wire. Further, 
adding the interaction effects has a considerable effect 
on this density difference, and now it is peaked both at 



X (nm) 
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FIG. 2: The difference between the single-particle densities 
of the many-body states in odd and even sector without (a) 
and with (b) Coulomb interaction. Majorana fermion proba¬ 
bility distribution for 71 of m without (c) and with (d) the 
Coulomb interaction. For all figures, we take fj, = 0.4 meV 
and B = 0.7 T. 


the center and at the ends of the wire, despite the fact 
that it was predicted 27| that the electron-electron inter¬ 
action will not have a considerable effect on the charge 
distribution due to the MFs. 

In Fig. [5] (c, d) the MF probability distribution is 
shown for of dill) without and with the Coulomb in¬ 
teraction taken into account. The figures for the MF 73 
of (ITT]) are inversion symmetric to the case of 7 ^ at the 
inversion point x = y = Q. It should be noted that the 
procedure of obtaining the MF wave functions outlined 
above is exact when in the topological superconducting 
phase the ground state is doubly degenerate, which as 
we saw above is not satisfied for the realistic system. 
Therefore the MFs which we obtain in our calculation 
do not generally satisfy the 7 ^ = 1 relation. This means 
that the MF wave functions obtained using the proce¬ 
dure above will only be approximately normalized for the 
system considered in this work. For the non-interacting 
case we get the normalization of the MF wave function 
to deviate from unity by about 4%, whereas for the in¬ 
teracting case it is around 10%. This difference between 
the non-interacting and interacting cases is related to the 
fact that for the interacting case the MF operator can 
have additional nonlinear terms in the expansion using 
the electron creation and annihilation operators and 
Cns [ij- As shown in Fig. [2] (c) for the non-interacting 
case the MF probability distribution is shifted toward the 
left part of the wire, although its maximum is close to 
the center of the wire than to the edge. Adding inter¬ 
action into the picture (Fig. [5] (d)) splits the high peak 
into two parts, one located at the left edge and the other 
located at the center of the wire. Therefore the MF wave 
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functions obtained with the procedure outlined in equa¬ 
tions (HUHl gives only approximate results in this exact 
diagonalization approach. 

In conclusion, we have considered the effect of long 
range Coulomb interaction on the topological supercon¬ 
ducting phase in a quasi ID semiconductor wire proxim¬ 
ity coupled to a s-wave superconductor using the exact 
diagonalization approach. By dealing with the problem 
using the exact diagonalization procedure we were able to 
treat the Coulomb interaction directly, without resorting 
to any approximation. We have found that adding the 
interaction into the system enlarges the region of the pa¬ 
rameter values, such as the magnetic field strength and 
the chemical potential, where topological superconduc¬ 
tivity can be observed. This is consistent with the previ¬ 
ous studies. In addition, we have also found that while for 


the magnetic field strength close to the phase transition 
point, the interaction lowers the bulk superconductivity 
gap, for moderate values of the magnetic field strength 
away from the transition point, the interaction actually 
enhances the bulk gap. This finding can be important 
for observation of the topological superconducting phase 
in the experiment due to the fact that the effect of the 
Coulomb interaction can be controlled by manipulating 
the charge density in the semiconductor wire. Finally, the 
present approach is quite general and can also be used 
to consider the properties of topological superconductiv¬ 
ity in a 2D semiconductor wire with long-range Coulomb 
interaction, and include additional features such as the 
orbital effects or the multichannel filling. 

The work has been supported by the Canada Research 
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